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Aims. We search for Polycyclic Aromatic Hydrocarbon (PAH) features towards young low-mass (T Tauri) stars and compare them with 
surveys of intermediate mass (Herbig Ae/Be) stars. The presence and strength of the PAH features are interpreted with disk radiative transfer 
models exploring the PAH feature dependence on the incident UV radiation, PAH abundance and disk parameters. 
Q ' Methods. Spitzer Space Telescope 5-35 (im spectra of 54 pre-main sequence stars with disks were obtained, consisting of 38 T Tauri, 7 Herbig 

\ Ae/Be and 9 stars with unknown spectral type. 
CZ2 . Results. Compact PAH emission is detected towards at least 8 sources of which 5 are Herbig Ae/Be stars. The 1 1.2 /jm PAH feature is detected 
in all of these sources, as is the 6.2 fim PAH feature for the 4 sources for which short wavelength data are available. However, the 7.7 and 8.6 
J> /im features appear strongly in only 1 of these 4 sources. Based on the 1 1.2 pm feature, PAH emission is observed towards at least 3 T Tauri 
• ^ ■ stars, with 14 tentative detections, resulting in a lower limit to the PAH detection rate of 8 %. The lowest mass source with PAH emission in 
|/\J our sample is T Cha with a spectral type G8. All 4 sources in our sample with evidence for dust holes in their inner disk show PAH emission, 
increasing the feature/continuum ratio. Typical 11.2 p.m line intensities are an order of magnitude lower than those observed for the more 
massive Herbig Ae/Be stars. Measured line fluxes indicate PAH abundances that are factors of 10-100 lower than standard interstellar values. 
Conversely, PAH features from disks exposed to stars with T eff < 4200 K without enhanced UV are predicted to be below the current detection 
limit, even for high PAH abundances. Disk modeling shows that the 6.2 and 11.2 pm features are the best PAH tracers for T Tauri stars, 
whereas the 7.7 and 8.6 /im bands have low feature over continuum ratios due to the strongly rising silicate emission. 
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1. Introduction 

Polycyclic Aromatic Hydrocarbons (PAHs) have been ob- 
served in a wide variety of sources in our own and external 
galaxies. Within the Milky Way PAHs are observed in the dif- 
fuse medium, dense molecular cloud s, circumstellar enve lopes, 
and (proto-)planetary nebulae (see IPeeters et alJ I2004I for a 
summary). A common characteristic of all of these sources 
is that they are exposed to copious ultraviolet (UV) photons. 
The UV radiation drives the molecules into excited electronic 
states, which subsequently decay to lower electronic states 
through a non-radiative process called internal conversion, fol- 
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lowed by vibrational emission in the available C-H and C-C 
stretching and bending vibrational modes at 3.3, 6.2, 7.7, 8.6, 
11.2, 12.8 and 16.4 fim. Thus, PAH molecules form an impor- 
tant diagnostic of UV radiation. 

In recent years, PAH emission has also been detected 
from disks around young stars i n ground-based and Infrared 
Space Observatory (I SO) spectra JVan Kerckhoven et al 1 2000 ; 
[ Honv et alJl200lMPeeters et alJl2002l Ivan Boekel et alj 12004 ; 
IPrzvgodda et alJ 120031 lAcke & van den AnckeJ l20ol . The 
PAH emission is thought to originate from the surface layer 
of a (flaring) d isk exposed to radiation from the central star 
(cf. models by iManske & Henninglll999t lHabartet alJl2004 
hereafter H04). Indeed, ground-based spatially resolved obser- 
vations show that the features come from regions with sizes 
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typical of that of a circumstellar disk (radiu s <12 AU at 3.3 
jum, <100 AU at 1 12 iim ) JGeers et all2005Uvan Boekel et alJ 
200-4 lHabart et alJl2005l) . Searches for PAHs in disks are im- 
portant because in addition to being a tracer of the strength of 
the stellar radiation field and disk geometry, the PAHs also af- 
fect the disk structure and chemistry. For example, th e high 
opacity of PAHs at FUV wavelengths (Mattioda et al. 2005) 
could significantly reduce the stellar UV radiation in the inner 
disk while photoionization of PAHs produces energetic elec- 
trons which are a major heating mechanism for the gas in the 
upper layers of the disk where the gas and dust temperatures are 
not w ell coupled ( Jonkheid et al. I l2004 iKamp & Dullemondi 
2004). 

Detections of PAH features also provide diagnostics of 
the presence of small grains in the surface layers of disks 
and the dust evolution through grain growth and dust settling. 
Evidence for grain growth has been found from modeling of 
the silicate features fro m Herbig Ae/Be (h ereafter HAe Be) 
Jvan Boekel et alJ 120041) and T T auri disks JPrzveodda et alJ 
2003t iKessler-Silacci et ail 120051 2006) and more indirectly 
from the mo deling of the observ ed H2 emission features from T 
Tauri disks (Bergi n et alJl2004l) . PAHs are considered to be on 
the small end of the size distribution of grains. An important 
question is whether PAHs have a different timescale for set- 
tling and/or growth compared to that of larger silicate/carbon 
dust grains. 

So far, most data obtained on PAHs refer to relatively bright 
features in the spectra of intermediate mass HAeBe stars. A 
recent ISO spectroscopic survey has detected PAH features, in 
particular the most frequently observe d 6.2 um feature, towards 
57 % of a sample of 46 H AeBe stars jAcke & van den AnckeJ 
l2004 . lMeeus et alJfcOOll) classified the ISO observed spectral 
energy distributions (SEDs) of these int ermediate mass young 
stars into two groups (group I and II). iDominik et alJ 12003) 
interpreted these SEDs in th e context of a passive d isk model 
with a puffed-up inner rim llDiillemonde talJ l200lh and pro- 
posed that group I sources with larger mid-infrared excess have 
flaring disks while group II s ources are consistent wit h small 
and/or self-shadowed disks. lAcke & van den AnckeJ (|2004) 
showed that the group I sources with strong mid-infrared rel- 
ative to near-infrared excess display significantly more PAH 
emission than the group II sources with weaker mid-infrared 
excesses, consistent with the idea that the PAH emission origi- 
nates mostly from the disk surface. 

The arrival of the Spitzer Space Telescope l|Werneretal 



20041) with the InfraRed Spectrograph (IRS) faouck et al 



2004) provides the opportunity to extend these studies to disks 



around fainter low mass T Tauri stars. For sources of spec- 
tral type G and later, the stellar UV field is orders of magni- 
tude weaker than for HAeBe stars, which will directly affect 
the PAH excitation and emission. On the other hand, enhanced 
UV radiation has been detected for some T Tauri stars (e.g., 
Costa et al. 2000, Bergin et al. 2003). Such an enhanced UV 
field should be directly reflected in the intensity of the PAH 
features if these molecules are present in normal abundances. 
Furthermore, ionized PAHs can be ex cited by less energeti c, 
optical photons from these cooler stars dMattioda et all2005l) . 



In this paper we present detections of PAH features toward 
T Tauri stars from our initial set of Spitzer IRS 5-38 /urn spec- 
troscopic observations that were taken as part of the Spitzer 
Legacy program "F rom Molecular Cores to Planet-Forming 
Disks" fevans et all2003l) (hereafter c2d). The c2d targets con- 
sist of a large number of sources with infrared excess in five 
of the nearest large star-forming regions: Chamaeleon, Lupus, 
Ophiuchus, Perseus and Serpens. In Paper I, the silicate 10 and 
20 fim features from pre-main sequence stars with disks are 
presented and analyzed (Kessler-Silacci et al. 2006). The data 
are used here to search for PAH features and address a num- 
ber of outstanding questions. For how many low mass stars can 
PAH features be found and how does this compare to HAeBe 
stars? Can limits be put on the abundance of PAHs and thus in- 
directly on that of the smallest grains? Which factors influence 
the appearance of PAH features in disks? Can we quantify any 
additional UV or optical radiation from the strength of the PAH 
features? 

Section[2]describes the sample selection, observations and 
reduction method. The results for the observed PAH fea- 
tures are presented in Sect.[3]and discussed in Sects. |4] and [5] 
Conclusions are presented in Sect.[6] 

2. Observations and data reduction 

Mid-infrared spectra were obtained for intermediate and low 
mass stars with circumstellar disks with the IRS aboard Spitzer 
as part of the c2d Legacy program. All targets were observed 
with the Short-High (SH) module (spectral resolving power 
of R ~ 600) covering the wavelength range 10-20 /mi and 
thereby potential PAH features at 1 1.2 and 12.8 fim. For a frac- 
tion of our sources, Short-Low (SL) observations (R ~ 100, 
5-14.5/im) were obtained as well, which cover the 6.2, 7.7, 
8.6, 11.2 and 12.8/mi PAH features. For the remainder of 
our sources, the SL observations are in the Guaranteed Time 
Observation (GTO) IRS program (J.R. Houck) and were at the 
time of paper submission not yet available. 

The or iginal sample sele ction of the c2d program is de- 
scribed in lEvans etafl J2003h . In our sample of sources ob- 
served to date there are 7 HAeBe stars, 38 T Tauri stars and 
9 sources without known spectral types. These include the 40 
T Tauri stars and 7 HAeBe stars in Kessler-Silacci et al. (2006), 
although we label 4 of their T Tauri sources here as "unclassi- 
fied" for lack of spectral type (RNO 15, IRAS 03446+3254, 
VSSG1, CK 4) and we label 1 of their HAeBe stars here as a 
T Tauri star since it has spectral type M6 (DL Cha). In addi- 
tion, we have added the c2d IRS spectra taken for 1 HAeBe 
star (HD101412), 1 T Tauri star, Sz 84, and 5 sources without 
known spectral types (EC 69, EC 88, EC 90, EC 92, RNO 91). 
Our current sample of 54 sources is therefore comparable in 
size to that of the HAeBe stars observed with ISO. We grouped 
the sampl e sources in in termediate and low mass young stars 
following lThe et alJ dl994l) . as HAeBe stars with spectral types 
of F7 and earlier, and T Tauri stars with spectral type F8 and 
later. A summary of the properties o f all disk sources can be 
found in lKessler-Silacci et alJ (I200(tL Paper I) and Merin et al. 
(in prep.). The sample includes the spectrum of an off-position, 
taken toward the Ophiuchus cloud at a position devoid of in- 
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frared sources, to compare with a typical interstellar medium 
PAH spectrum. 

All spectra were extracted from the SSC pipeline ver- 
sion S 12.0.2 BCD images, using the c2d reduction pipeline 
(Kessler-Silacci et al. 2006 and Lahuis, in preparation). The 
processing includes bad-pixel correction through interpolation 
using a source profile fit in the cross-dispersion direction. The 
source profile fitting also gives an estimate of the local sky con- 
tribution in the high-resolution spectra. The extracted spectra 
are defringed using the IRSFRINGE package developed by the 
c2d team and individual orders are in some cases corrected by 
small scaling corrections (<5 %) to match the order with the 
shortest wavelength. Sky subtraction is applied and in a few 
cases the SH module is scaled down (< 5 %) in flux to match 
the SL spectrum. 

3. TTauri stars with PAH features 

3.1. Identification of PAH features 

Among the 54 disk sources observed to date, a relatively small 
sample of 8 sources shows one or more clear emission feature 
that we attribute to PAHs (Fig.^. Table[2summarizes the char- 
acteristics of these sources. They consist of 3 TTauri stars and 
5 HAeBe stars. Besides the PAH features, some sources show 
the H2 S(2) line, whereas the [N e II] 12.8 fim line is detected 
toward T Cha. Unlabelled narrow features are likely spurious. 

Disk sources with PAH emission features were identified 
as follows. When both SL and SH spectra are available, we re- 
quired detection of both a clear 1 1 .2 fim feature and a clear 6.2, 
7.7 and/or 8.6 /im feature. For sources where only a SH spec- 
trum is available (A > 10/mi), the identification is more criti- 
cal and largely relies on both the presence of a clear feature at 
11.2 fim and its shape. We compared potential 11.2 fim features 
with that observed in the off-position spectrum and selected the 
sources with the lowest residuals, see Sect. 14. 31 

This method can introduce a bias against sources with 
mixed crystalline silicate and PAH features, because crystalline 
forsterite has a characteristic feature at almost the same cen- 
tral wavelength as the 1 1 .2 fim PAH feature. To distinguish be- 
tween these two possible assignments, we searched for other 
crystalline silicate features at wavelengths longer than 11.2 fim 
(e.g. 16.2, 18.9, 23.7 and 33.6^/m). The large wavelength cov- 
erage out to 35 fim is a significant advantage compared with 
ground-based data. Three of the five HAeBe and 2 of the 3 T 
Tauri sources show (tentative) evidence for crystalline silicates 
at longer wavelengths, indicated in Table fj] (Kessler-Silacci et 
al. 2006). In these sources, a contribution from crystalline sili- 
cates to the observed 1 1 .2-1 1 .3 fim feature cannot be excluded. 
The sources illustrate well the difficulties inherent in the iden- 
tification of PAH emission features at 11.2 fim when lacking 
any information on the presence/absence of other characteris- 
tic PAH features at shorter wavelengths. We note also that the 
11.2 fim PAH feature can be blended with the broad amorphous 
silicate feature whose strength and spectral width vary with 
grain size. This may have led us to miss sources with weak 
11. 2 fim emission from PAHs (see discussion in Sect. l3.3l l. al- 
though the characteristic shape of the 11.2 fim PAH feature can 
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Fig. 1. Spitzer IRS spectra of sources with PAH features, com- 
prised of the SL (5-10 fim) and SH (10-20 fim) modules. The 
location of PAH features is indicated with markers at 6.2, 7.7, 
8.6, 11.2 and U.Sfim. 



help in distinguishing between PAHs and silicates, as discussed 
in Sect.l4~3l 



Finally, sources with PAH features but showing silicate 
and/or ice features in absorption have been excluded from the 
sample presented here. For comparison, Fig.|2]includes obser- 
vations of the off-source spectrum and of two more late type 
sources, Haro 1-17 and VSSG1 which were selected on the 
above mentioned criteria but for which background extraction 
shows that the PAH emission features are fully associated with 
the extended cloud emission (see below). 
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Table 1. Sources with potential PAH emission and their characteristics 
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a SL = SL1 + SL2 

b average from 3 measurements 

c average from 19 measurements 

d PAH feature detected in background spectrum, not associated with the source 
c label 'e' added to spectral type here, based on Ho- > 0, not taken from reference. 

f distance uncertain m m m m m 

References for distance; s pectra l type; Hor: l: lvan den Ancker et alj <1998t). 2:lAlcala et alj <1993l). 3 : Enoch et al. 2005. 4: | Fernandez et alj 
19951) . 5: ICohen & KuhlNl979l). 6: As sumed distance to O ph cloud jde Geus et alll989h . 7- lPrato et aljEool) . 8:lAcke & van den AnckeJ 
120041). 9: iDunkin et a ll <1997h. 10: lAcke et alJ J2005h. 11: Assumed dist ance to Taurus-Auriga cloud iKenyon et alJ <1994li . 12: 
Hernandez et alj j2004). 13: iMora et alj |200ll). 14: IStraizvs et alj <1996h . 15: iFinkenzeller & Mundtl Jl984l) . 16: iThe et alj 7l994l) . 17: 
Hould <1978p 8: Ulcala et alJJl995h . Lg JlVlartin et alJll998l) . " 

Note: sources in the bottom portion of the Table have spectra with PAH features that are fully attributed to background emission. 



3.2. Spatial extent of PAH emission 

An important question is whether the PAH emission originates 
primarily from the star+disk system or from an extended neb- 
ulosity around the star. To properly answer this question, either 
spatially resolved spectra or images in specific PAH band filters 
(both at feature wavelength and slightly off-peak for determi- 
nation of strength and extent of continuum emission) with suf- 
ficient (subarcsec) spatial resolution are required. Since such 
data are lacking for our sample we use the extracted back- 
ground and long-slit spectra to provide constraints on the extent 
of the emission. 

The SL spectra are taken using long-slit spectroscopy and 
the 2D spectral images can be used to determine if the features 
are seen extended along the entire width of the slit. An exam- 
ple 2D spectral image of module SL1 is shown for RR Tau in 
Fig.|3j where no extended emission is seen along the slit. The 
pixel size of the SL module is 1.8", which means that any fea- 
ture originating from a region smaller than 2 pixels (e.g. 540 
AU at 150pc) will be spatially unresolved. A similar lack of 
extended emission in SL1 is seen for all of our PAH sources 
with the exception of Haro 1-17. Thus, the PAH emission for 
these sources is constrained to originate from a region of at 
most 3.6". 

The SL spectra presented have been corrected for the esti- 
mated background contribution (see Sect.|2j. In 4 of the 5 SL 
spectra (LkHa 330, RR Tau, HD 101412 and VV Ser), the fea- 
tures remain after background subtraction and are concluded to 
be associated with the source. However, for 1 source, Haro 1- 



17, the sky-corrected SL spectrum shows no PAH features (see 
Fig. |3J . The PAH emission seen towards this source is con- 
cluded to be entirely due to background emission. 

A number of the sources were only observed in SH for 
which the width of the slit is too small to directly extract 
the continuum flux outside the source profile. Here, the back- 
ground contribution is estimated from fitting a standard star 
source profile plus a flat continuum flux to the measured source 
profile. This background subtraction is more difficult, since the 
small slit length of only 1 1 .2" (5 pixels) makes the source pro- 
file fits less accurate. 

For all presented sources, the SH source and sky spectra 
were inspected for PAH features. Two SH background spectra, 
one with extended PAH features (VSSG1) and one without (RR 
Tau) are shown in Fig.|4] The extracted background spectra for 
SH show no PAH features for 8 of the 10 sources in our sam- 
ple with PAH emission. The two sources with features in the 
background spectra are Haro 1-17 and VSSG1. In both cases, 
removing the background emission entirely removes the 11.2 
/i/m PAH feature from the source spectrum (Fig. |2j. It is con- 
cluded that for these two sources the PAH features are fully due 
to background emission. 

The environment around one source, the Herbig Ae star 
VV Ser, is discussed extensively by Pontoppidan et al. (2006). 
Based on IRAC and MIPS images, it is shown that VV Ser is 
surrounded by a bright and extremely large (~ 6') nebulosity 
emitting at 8, 24 and 70 /vm, which is not seen in near infrared 
images. They conclude that the emission is due to a mix of 
quantum-heated PAHs and Very Small Grains (VSGs) present 
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Fig. 2. Spitzer IRS spectra (black line) of sources with PAH 
features that are associated with extented cloud emission, com- 
prised of the SL (5-10 /mi) and SH (10-20 /mi) modules. The 
extended emission corrected spectrum is shown in grey. The 
location of PAH features is indicated with markers at 6.2, 7.7, 
8.6, 11.2 and 12.8 /mi. 





Fig. 3. 2D spectral image of modules SL2 (left panel) and SL1 
(right panel) for RR Tau. The quantity shown is the ADU in 
units of e~s _1 . The positions of the 6.2, 7.7, 8.6 and 11.2 /mi 
PAH features are indicated with labels. The horizontal bar at 
bottom right of the left panel and top left of the right panel 
indicates a spatial extent of 20". 



in the low-density nebula. The PAH emission lines in the spec- 
trum must originate from within 1 .8" - 2.35" or 450 - 600 AU 
(half-slit width SL-SH). In their best-fit model they require a 
central cavity of 15000 AU, so the extended nebulosity is not 
expected to account for the observed PAH features although 
they note that it cannot be ruled out that a small clump of PAH 
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Fig. 4. IRS SH spectra for VSSG 1 and RR Tau. Bottom 
(black): extracted sky spectrum. Middle (grey): source spectra 
corrected for background. Top (black): source spectra. The sky 
spectra are shifted by +0.5 and + 1 .0 Jy respectively for purpose 
of clarity. 

material is present nearer to the star. Here we assume the emis- 
sion is from PAHs in the circumstellar disk. 

In summary, we conclude that the PAH emission from most 
of our sources does not originate from extended diffuse fore- 
or background emission and must instead originate from the 
observed young stars with disks. 

3.3. Statistics 

Within our current c2d sample, clear PAH features are detected 
in most (5 out of 7) HAeBe stars while only 3 out of 38 T 
Tauri stars show features consistent with the presence of PAH 
molecules. Interestingly, the PAH detection rate is 100% for 
the 4 sources (HD135344, SR 21N, LkHor 330 and T Cha) with 
SEDs characteristic of cold disks, i.e., sources with SEDs that 
lack excess emission in the 3-13 /mi region, indicative of an in- 
ner hole in the dust disk (Brown et al., in prep). These 4 sources 
are also the 4 lowest mass sources (spectral type F4 - G8) with 
PAH detections. 

The c2d sample of PAH detections is biased towards 
sources with either a strong 1 1 .2 /mi feature or multiple PAH 
features from SH and SL observations. This excludes potential 
sources with weak PAH features for which only SH observa- 
tions are available now: since we cannot assign the origin of 
this 11.2 /mi feature to either PAH or crystalline forsterite, these 
sources are, for now, excluded. This includes 17 sources with 
tentative 1 1.2 yum detections: 14 T Tauri stars, 1 Herbig Ae star 
and 2 unclassified sources. The 14 T Tauri stars are DoAr 24E, 
EC 82, GW Lup, GY 23, HT Lup, Krautter's Star, RU Lup, SX 
Cha, SY Cha, SZ 73, VW Cha, VZ Cha, V710 Tau (binary) and 
WX Cha. Future SL data will be able to confirm or dismiss the 
presence of PAHs in these sources. 

Finally, in the current sample there are 28 sources with 
no clear 11.2 /mi PAH feature, consisting of 1 Herbig Ae star, 
21 T Tauri stars and 6 sources with no known spectral type. 
From our present sample, the lower limit to the detection rate 
of PAH emission features toward T Tauri stars is about 8%. 
This detection rate goes up to 45 % if tentative detections are 
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Fig. 5. Left panel: Blow-up of the SH spectra around the 
continuum-subtracted 11.2 pm PAH feature, normalized to the 
fitted peak flux (black line). Overplotted in light grey for com- 
parison is the SH spectrum of the off-position. Right panel: 
Plot of the difference between the source and the off-position 
spectra. 



included. PAH emission is only detected toward the 3 more 
massive T Tauri stars of spectral type G of our sample, which 
consists mostly of K-M type stars. 

The rather small fraction (8 %) of low mass stars in our 
c2d sample with PAH emission features, based on the 1 1.2 pm 
P AH feature, contrasts with the la rge fraction (57 % according 
to lAcke & van den Anckerl 12004) for intermediate mass stars 
detected with ISO based on the 6.2 pm PAH feature. If only the 
11.2 pm features are considered, their detection rate drops to 
48 % based on their Table 3. However, this potentially includes 
sources with only crystalline silicate emission. 




5.5 6.0 6.5 7.0 5.5 6.0 6.5 7.0 

Wavelength (pm) 



Fig. 6. Left: Blow-up of the Spitzer-IRS low resolution spectra 
around the 6.2 //m PAH feature. A simple fit of the continuum 
flux is plotted with a dotted line. Right: Continuum subtracted 
spectra, normalised to the peak flux of the PAH feature. 

4. Analysis of PAH features 

4.1. Overview of detected PAH features 

Table |2 summarizes the detected PAH features and measured 
line fluxes. The spectrum of RR Tau nicely shows all the main 
PAH bands detectable in the IRS spectral window at 6.2, 7.7, 
8.6, 11.2, 12.8 pm (Fig.^ and even 16.4pm (not shown). In 
contrast with the other sources, it has clear 7.7 and 8.6 pm fea- 
tures, attributed to C-C stretching transitions and C-H in plane 
bending transitions. These features are either less obviously de- 
tected above the silicate continuum or simply absent in the 4 
other sources with SL spectra. LkHo' 330 is a clear example 
of a spectrum with a 6.2 /mi C-C stretching feature and a 1 1.2 
pm C-H out-of-plane bending feature but no 7.7, 8.6 nor 12.8 
pm features. The 7.7 and 8.6 pm features are ge nerally f ound 
to be well correlated with the 6.2 pm feature JPeeters et alJ 
2002) and the ir absence is th us puzzling in our high sensitivity 
spectra. ISooon et"aD ((2002) show that silicate absorption can 
strongly mask the 8.6 pm PAH feature ; in our case the silicate 
is in e mission, however. According to lAcke & van den AnckeJ 
(2004) (see their Table 3), four HAeBe stars observed with ISO 
similarly show 6.2 and 1 1.2 pm emission but no features at 7.7 
or 8.6 pm. Two of these sources, HD 163296 and VV Ser, are 
also in the c2d sample (the other two are HD 142666 and HD 
144432). This absence is further discussed in Sect. 15.51 in the 
context of a disk model which demonstrates the lower contrast 
of the 7.7 and 8.6 pm features with respect to the continuum 
emission from the disk. 

4.2. Line flux determination 

To determine the strength of the PAH features, the continuum 
emission needs to be subtracted. As a simple approximation 
we derive a local pseudo-continuum by fitting a 2D polyno- 
mial to the spectrum around the individual PAH features, where 
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Line flux feature/continuum 



Name 


6.2/jm 
(SL) 


11.2yum 
(SL) 


11.2pm 
(SH) 


12.8pm 
(SH) 


6.2 pm 
(SL) 


1 1.2 pm 
(SH) 


cryst. 
sil.° 


LkHo- 330 


1.7 x 10~ 15 


7.7 x 10~ 16 


5.1 x 10- 16 


< 2.5 x 10 16 


1.21 


1.12 


T 


RR-Tau 


9.6 x 10~ 15 


5.0 x 10~ 15 


4.5 x 10~ 15 


1.3 x 10~ 15 


1.71 


1.57 


Y 


HD 98922 




a 


1.4 x 10~ 14 


4.0 x 10~ 15 




1.06 


Y 


HD 101412 


4.0 x 10- 15 


2.2 x 10- 15 


1.5 x 10~ 15 


3.0 x 10~ 16 


1.19 


1.08 


yd 


TCha 




a 


3.3 x 10~ 16 


1.1 x 10~ 16 




1.15 


N 


HD 135344 




a 


1.2 x 10~ 15 


1.1 x 10~ 16 




1.19 


N 


EM* SR21 N 




_a 


4.0 x 10- 15 


4.1 x 10~ 16 




1.32 


T 


VV Ser 


2.6 x 10~ 15 


3.1 x 10~ 15 


2.3 x 10~ 15 


1.2 x 10" 15 


1.07 


1.07 


Y 


Off-position 1 
VSSG 1 b 
Haro 1-17 b 


a 
_a 

4.5 x 10~ 16 


a 
_a 

< 2.5 x lO" 16 


1.4 x 10~ 16 
8.6 x 10- 16 
< 2.5 x 10-' 6 


5.5 x 10- 17 
2.4 x 10~ 16 
< 2.5 x 10~' 6 






Y 
Y 



a no SL spectra available 

b PAH feature detected in background spectrum, not associated with the source 

c crystalline silicates detected in either 28-29, or 33-35 /an, from Kessler-Silacci et al. 2006, Table 2; "Y" if detected, "N" if not detected, 
or "T" if the identification is tentative. 
d derived in this study 



the continuum is selected by hand. For the 1 1 .2 pm feature, 
a polynomial is fitted to the emission at 10.5-11.0 and 11.8- 
12.2yum. The continuum emission below the 6.2 pm emission 
is estimated between 5.5 and 7.1 pm. We do not include the 7.7 
and 8.6 pm features in Table[2]because these do not appear sig- 
nificantly in our sources with SL data, with the exception of 
RR Tau. The 12.8 pm line fluxes extracted from SL are within 
the uncertainty consistent with those extracted from SH. 

The continuum-subtracted features are integrated between 
fixed wavelengths, in particular between 6.0 and 6.6 /mi for the 
6.2 pm PAH feature, between 10.9 and 1 1 .6 /im for the 1 1 .2 pm 
feature (see Fig. IB.ll in the online material) and between 11.8 
and 13.2pm for the 12.8pm feature. The resulting continuum 
subtracted features are shown in Figs.l5ll6land lB.2l 

The measured line fluxes are summarized in Table [2] 
Thanks to the increased sensitivity of Spitzer, our derived line 
fluxes for clearly detected PAH features are an order of mag- 
nitude lower that what was previously possible with ISO. Our 
weakest detected feature is the 11.2 pm PAH feature in T Cha 
with a line flux of 3.3 x 10~ 16 W trT 2 . A mean 3cr sensitivity 
limit of 2.5 x 10~ 16 W rrT 2 is derived from the noise determi- 
nation in the continuum adjacent to the PAH features, though 
this limit varies somewhat from source to source, depending on 
differences in the S/N of the reduced spectra and on the pres- 
ence of residual reduction artifacts for a few cases. For a small 
number of sources the sensitivity limit reaches a few x 10~ 17 
WnT 2 . 

Three HAeBe stars in our sample — HD 135344, RR Tau 
and VV Ser — have previousl y been observed with ISO, albei t 
with much lower S/N ratio lAcke & van den AnckeJ 12004). 
For RR Tau, our derived line flux for the 6.2 pm feature agrees 
within ~ 5 % with the ISOPHOT-S spectra, while the 11.2 pm 
feature is larger by about a factor 2 in the IRS spectrum. For 
HD 135344, we derive a slightly higher 11.2pm line flux than 
the ISO upper limit. For VV Ser, our derived 6.2pm line flux 



is about 4 times weaker, whereas our 11.2 pm detection is con- 
sistent with the ISO data. 

Sloa n et all J2005I) have presented Spitzer SL observations 
for 4 HAeBe stars with spectra showing PAH features but no 
silicate dust features, among which HD 135344 is included in 
our sample. They report clear 6.2, "7.9", 11.3 and 12.7 pm 
features for all of their sources (HD 34282, HD 135344, HD 
141569, HD 169142). For HD135344 the 7.9 pm feature is 
weaker and broader and the 8.6 pm feature is absent. Their de- 
rived line flux for the 11.2 pm PAH feature is consistent with 
ours within 10%. 

4.3. Comparison of PAH features 

In a previous study of a diverse sample of interstellar and 
circumstellar sources, plan etary nebulae, reflection ne bulae 
HII regions and galaxies Jvan Diedenhoven et alJ Eo04l) . the 
1 1 .2 pm PAH feature of almost all YSO's, non-isolated HAeBe 
stars and HII regions was found to have a similar asymmet- 
ric profile with a FWHM of ~ 0. 17 pm, and a peak wavelength 
in the range of 1 1 .2 - 1 1 .24 pm. The single isolated HAeBe 
source in their sample, HD 179218, shows a broader 11.2pm 
PAH feature with a peak wavelength of ~ 1 1 .25 pm. Our peak 
position lies in the range 1 1 .25 - 1 1 .32 pm + 0.03 pm, but this 
determination is influenced by the uncertainty in the contin- 
uum fit, and the shift to longer wavelengths may be partially 
explained by the presence of a 11.3 pm feature from crystalline 
silicates. 

Figure |3 shows a comparison of the continuum-subtracted 
11.2 pm features from all SH spectra. The off-source PAH fea- 
ture, which is clearly not contaminated by silicate emission, is 
included. Both the source and the off-source features are nor- 
malized to the peak flux. In the right panel of Fig.[3]the differ- 
ence between the two features is shown. 
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For T Cha, LkHa 330, SR 21 N and HD 135344, the 
11. 2 /mi feature is broader than, as well as redshifted with 
respect to, the off-position feature. Of these, T Cha and HD 
135344 show no evidence for crystalline silicates at 28-29 and 
33.6 yum (Kessler-Silacci et al. 2006), while LkHar 330 and SR 
21 N show tentative crystalline features. Thus, the presence 
of a 11.3 fim crystalline silicate feature cannot readily explain 
the broadening of the measured feature for these sources. Such 
a broad shape has been seen before in the planetary nebula e 
IRAS 17047-565 and IRAS 21282+5050 felony et alJl200ll) . 
IPechet al.N2002l) proposed anharmonicity as an explanation of 
the broadening and used a PAH emission model to fit the 1 1 .2 
fim feature of IRAS 21282+5050 with a combination of the 
fundamental (v = 1 — > 0) and hot bands (v = 2 — > 1 and 
v = 3 — > 2) of the transition. These hot bands would point to 
very hot PAHs being present, presumably in the innermost part 
of the disk where the radiation field is strongest. 

For the HAeBe sources RR Tau, VV Ser, HD 98922 and 
HD 101412, the peak wavelength of the 11.2 PAH feature is 
very similar to that of the off-position feature. For RR Tau, the 
PAH feature compares very well with the off-source feature, 
both in shape and peak position, showing little to no residual 
after subtraction. For the other 3 sources, subtracting the off- 
position feature leaves several residual features, hinting at the 
presence of crystalline features, which are also seen at longer 
wavelengths. 

5. PAH emission from disks 

5.1. Disk model 

The strength of the PAH emission features is known to depend 
on the strength of the UV and optical radiation field, but in 
disks several additional parameters can affect the appearance of 
the PAH features. Here we address the question as to how the 
PAH features are affected by the spectrum of the central source, 
the PAH abundance and the flaring geometry of the disk. A 
related question is why no PAH features are seen in at least 
half of our sample of T Tauri disks, in contrast to the findings 
for HAeBe stars. 

We use the 3- dimensional Monte Carlo rad iative trans- 
fer code RADMC JPullemond & Dominid EooH) . for which 
a module to treat the emission from quantum-heated PAH 
molecules and Very Small Grains has been included. This mod- 
ule will be described in detail in Dullemond et al. (in prep.), 
but a rough description has been given by Pontoppidan et 
al. (2006). A template model is set up using the following 
model parameters. A Kurucz model spectrum is taken for the 
central star with T e s = 10000 K, but the stellar parameters (lu- 
min osity, stellar rad ius) were chosen from evolutionary tracks 
bv lSiess et alJ d2000h for an age of 3 Myr. The disk is modeled 
with Mdisk = 1 X 1O~ 2 M , an inner radius set by a dust evapo- 
ration temperature of 1300 K, and an outer radius of 300 AU. 
The disk is flaring, with the vertical pressure scale height (in 
units of radius) at the inner rim H p IR m = 0.02 and at the outer 
edge H p /R out = 0.14. The disk is modeled to be close to face-on 
(i = 5.7°), to maximize the strength of the PAH features. The 
spectra are scaled to a distance to the observer of 150 pc. 



The PAH emission is calculated for an equal mix of 
neutral and singly ionized C100H24 molecules, adopting the 
lDrainefeLil f200n PAH emission model, using the "thermal 
continuous" approximation. Multi-photon events are included 
for the PAH excitation, f ollowing the method outlined by 
ISiebenmorgen et all Jl992h . Enhancement factors for the inte- 
grated cross sections of the 6.2, 7.7 and 8.6 um band s (Ef, ? = 3, 
£7.7 = 2, £8.6 = 2) as suggested by Li&Draine (2001) are 
taken into account, as also implemented in H04. The actual ion- 
ization state of PAH, which varies as a function of location in 
the disk, can affect the 7.7 and 8.6 /im features relative to 1 1.2 
and 3.3 fim but this will be explored further in future mod- 
els which include a full ionization balance of multiple PAH 
species i n disks (Visser et al. i n prep). We include the opaci- 
ties from Mat tioda et alJ d2005h for near-infrared wavelengths. 
Several model tests were performed, which are presented in 
AppendixlAl including a comparison with H04. 

PAH destruction is expected to occur in a strong UV ra- 
diation field, when the PAH molecules, through multi-photon 
events, absorb mor e than 21 eV in an interva l shorter than their 
cooling timescale ( Guhathakurta & Draine 1989). If the PAH 
destruction happens on a shorter timescale than the lifetime of 
the disk, this can have an effect on the PAH abundance in, and 
expected PAH emission from, the inner region of the disk. This 
effect will be stronger for smaller (N c < 50) PAHs, and de- 
pends on the assumed temperature at which the PAHs disso- 
ciate. Model calculations by Visser et al. (in prep.) show that 
for PAHs with N c = 100 the lifetime is larger than the life- 
time of the disk, and therefore the PAHs are kept at constant 
abundance throughout the entire the disk. Tests with remov- 
ing PAHs inwards from a destruction radius ^pAH.in (by setting 
their abundance to zero) show that the results do not depend 
sensitively on the choice of /?pAH,in when it is of order 1 AU. 

In our template model, we use the PAH abundance of H04, 
who adopt a fraction of 23 % of available carbon being locked 
up in C100H24 PAHs, corresponding to a carbon abundance with 
respect to hydrogen of 5 x 10 . For a single PAH species with 
Nc = 100 (C100H24) this leads to an abundance of 5 x 10~ 7 . 
In our model the abundance is set as a fraction of dust mass. 
Assuming C100H24 (molecular mass: 2.04 x 10~ 21 g) and a dust- 
to-gas ratio of 1 to 100, the PAH abundance of 5 x 10~ 7 with 
respect to hydrogen corresponds to a mass fraction of 0.061 
gram of PAH per gram of dust, of which 50% ionized and 50% 
neutral. Such an abundance is at the high end of that inferred for 
gener al interstellar clouds dCesarskv et aljEoOOl lHabart et alJ 
2001). Indeed, the parameters chosen in our model and that of 
H04 maximize the PAH emission. 

5.2. Dependence on spectral type 

Since UV and optical radiation incident on the disk surface pro- 
vides the main PAH excitation mechanism, the PAH features 
must depend on the spectral type of the illuminating star. To ad- 
dress this influence our template disk model is calculated for a 
range of Kurucz stellar spectra, with 7/ eff = 10000, 8000, 6000, 
5000 and 4000 K, corresp onding to spectral types AO, A6, GO, 
K2 and K7, respectively jGrav & Corballvlll994 . Stellar pa- 
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Fig. 7. Top panel: Model SEDs of a star+disk with PAHs, for 
different spectral types of the star, with T e ff = 10000, 6000, 
5000 and 4000 K at a distance of 150 pc. Middle panel: PAH 
line flux of the various features for the corresponding models, 
including also r e jf = 8000 K. Bottom panel: Feature to con- 
tinuum ratio of the PAH features. The dotted line indicates our 
3cr observational limits. 



ramet ers are again taken from evolutionary tracks JSiess et all 
2000) for an age of 3 Myr. The disk parameters were mod- 
ified slightly in 1 iteration to ensure hydrostatic equilibrium 
throughout the disk, except for the inner rim. The model SEDs 
for the different central stars and absolute continuum subtracted 
fluxes for the PAH features are presented in Fig.0 A blowup 
of the standard model spectrum is shown in Fig. [8] 

Since the absolute strength of the PAH features scales fore- 
most with the total radiation that is absorbed by the PAHs, it 
also depends on disk parameters that are unrelated to the PAHs 
themselves. To evaluate the role of the disk continuum, we 
present feature/continuum ratios in Fig.Qas well. 

It is seen that with decreasing T e fi both the line flux as well 
as the feature/continuum ratio decrease for the 3.3, 6.2, 7.7, 8.6 
and 11.2/im features (Fig.Q. The rate of decrease in line flux 
is similar for all features, as is that of the feature/continuum ra- 




5 10 15 

Wavelength (yum) 

Fig. 8. Blowup of the standard model spectrum of a star+disk 
with PAHs, for T e g = 6000 K with template parameters (dark 
solid line), without PAHs (dash-dotted line), with only neutral 
C100H24 PAH (light solid line). Note that the 11.2 and 6.2 //m 
features are clearly visible on top of the continuum, whereas 
the 7.7 and 8.6 //m features are masked by the rising 10 yum 
silicate feature. 



tio except for the 3.3 yi/m feature. The latter feature disappears 
more rapidly with decreasing r ef f because the continuum emis- 
sion at 3.3 yum decreases at a slower rate compared to the longer 
wavelengths. For all T e s considered, the 8.6 yum feature has the 
lowest feature to continuum ratio, from < 7 to < 3% , which 
can be clearly seen in Fig. [8] When the PAH molecules are as- 
sumed to be 100% neutral C100H24, the feature/continuum ratio 
of the 6.2, 7.7 and 8.6 yum bands decreases by a factor ~2-3, 
while that of the 3.3 and 1 1.2 //m features increases by a factor 
-1.5-3 (Fig. U. 

These models show that the 6.2 and 11.2 fim features are 
the most suitable tracers for the presence of PAHs in disks 
around late type T Tauri stars owing to their relatively high 
feature/continuum ratio, and this qualitatively agrees with the 
fact that these two features are always present in our T Tauri 
stars with PAH detection. Quantitatively, Fig.0shows that, for 
the current model assumptions, the PAH features become dif- 
ficult to observe for T e ff < 4200 K (later than K6) even with a 
relatively high PAH abundance of ~ 6% of the dust mass in the 
entire disk and large flaring angles. 

The observed 11.2 ^m line fluxes have all been rescaled to 
a distance of d — 150 pc, and are plotted in Fig.|5]together with 
the models. Also included is our typical 3<r limit for a source 
at 150 pc. For comparison w e added a small subsample of 1 2 
Herbig Ae stars, selected from Ac ke & van den A ncker 
and rescaled to 150 pc. Stars with unknown distance and/or 
poorly known spectral type were excluded. The observed fea- 
ture/continuum ratios for the two strongest features, 6.2 and 
11.2 fj.m, are compared with models in Fig. 1101 

Figure[9]shows that our template model (solid line) predicts 
a larger /pah than observed, with the exception of one source 
(HD 98922) for which the distance is uncertain. A lower ob- 
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Fig. 9. Strength of 11.2 ym PAH feature 7 PA h (scaled to a dis- 
tance of 150 pc) versus T e ff. Black solid line represents our 
template model, based on Fig. [71 dashed and dash-dotted lines 
represent our models with a lOx and lOOx lower PAH abun- 
dance respectively. Black '*' symbols are for c2d sources. Grey 
diamonds are ISO upper limits of the 11.2 yum feature strength 
for Herbig Ae sources, grey '+' symbols are ISO detections. 
c2d sources are labelled as follows. 1: TCha; 2: LkHa330; 3: 
SR21 N; 4: HD 135344; 5: RRTau; 6: VVSer; 7: HD 101412; 
8: HD 98922. The dotted grey line indicates our typical 3cr sen- 
sitivity limit of 2.5 x 10~ 16 W irT 2 for sources at d — 150 pc. 
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served 11.2 yum flux can be caused by a number of conditions 
such as a lower PAH abundance, a smaller flaring angle or a 
disk orientation close to edge-on (see below). Included in Fig. [9] 
are two additional model runs, where the PAH abundance has 
been lowered by a factor of 10 (dashed line) and 100 (dash- 
dotted line). Most of the observations fall within the predictions 
of these two models. Only the model with a lOOx lower PAH 
abundance predicts a 1 1 .2 yum feature strength for G-type stars 
below our Spitzer detection limit at d = 150 pc. One source, 
T Cha, is detected despite the fact it falls below this formal 
detection limit because it has a distance of d — 66 pc. 

Figure Qo] indicates that the template model fits the ob- 
served 11.2 yum feature/continuum ratios for the 4 low mass 
sources, even though it overpredicts /pah- The model with 
lOx lower PAH abundance fits the feature/continuum ratio of 
three of our HAeBe stars. However both models with lOx and 
lOOx lower (below plot limit) PAH abundance predict a fea- 
ture/continuum ratio below our Spitzer detection limit for spec- 
tral types below GO. This suggests that the small number of T 
Tauri detections (3 out of 38) are presumably outliers, with ab- 
normally high feature/continuum PAH features. This discrep- 
ancy will be further discussed in Sect. 15.51 and is most likely 
due to an abnormally low continuum in these sources. 

5.3. Additional UV radiation and relation with Ha 

The models presented in Sect. l5.2l include only the stellar radi- 
ation, but not any additional sources of UV. Excess UV radia- 
tion compared with the stellar photosphere has been observed 
from at least some T Tauri stars (e.g., Herbig & Goodrich 1986, 
Costa et al. 2000, Bergin et al. 2003). From a constructed com- 
posite FUV spectrum, derived f rom two K7 sources represen- 
tative of low-mass T Tauri stars, Ber gin etal J 120031 and refer- 
ences therein) find an overall FUV continuum flux at r — 100 
AU on the order of a few hundred times the Habing field. The 
effect of this additional UV is to shift absolute feature strengths 
to those appropriate for higher T e ff, close to T e ff ~ 10000 K. If 
a large fraction of T Tauri stars would have such additional 
UV radiation this leads to the question why not more sources 
were detected, since the feature strength should be more than 
sufficient if PAHs are present at our template abundance. This 
would further strengthen the evidence for low PAH abundances 
in a large fraction of T Tauri disks. 

The origin of excess UV radiation from T Tauri stars is not 
fully clear, but is usually thought to originate from the shock 
in the magnetospheric accretion column associated with ma- 
terial falling from the inner disk onto the stellar surface. One 
indication of accretion activity is the strength of t he Ha emis- 
sion l ine, measured in Ha equivalent width (EW) dCabrit et alJ 
1199(1 Literature values of Ha EW measurements are listed in 
Table [2 In our sample, LkHa 330 is a classical T Tauri star 
(Ha EW > 10 A), while T Cha should be classified as a weak 
line T Tauri star (Ha EW < 10 A). We find no correlation be- 
tween the 1 1 .2 /im PAH feature strength and Ha EW. However, 
single measurements of Ha EW may not necessarily be a reli- 
able tracer of the average accretion lumi nosity. As an examp le, 
for T Cha the quoted Ha EW of 2 A bv lAlcala et ail ill 9931) is 
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indicative of a low accretion rate, while ISvlvester et all 0996) 
have later classified it as a YY Orionis star with strong UV con- 
tinuum emission and substantial variability on short timescales. 
Another problem is that not all of this UV radiation from the 
accretion shock may reach the disk, but can be absorbed by the 
over lying accretion column and/or any bi-polar outflow mate- 
rial (Alexande r et alJ2005l) . 

An alternative interpretation of the enhanced UV radia- 
tion is stellar activity, especially at later evolutionary stages 
when the accretion becomes less importan t. This mechanism 
has been used by Kamp & Sammar (2004) to model the UV 
field around young G-type stars. 

5.4. PAH abundance 



The discussion in Sect. l5.1l has illustrated some of the effects of 
the PAH abundance on the PAH features. As seen in Figs.[9]and 
1101 the PAH features decrease in line flux and feature to contin- 
uum ratio with decreasing PAH abundance, with a similar trend 
for all the PAH features. Decreasing the PAH abundance also 
affects the overall SED, increasing the continuum radiation at 
wavelengths longer than 2 fim and at UV wavelengths around 
0.2yum. 

For sources with detected 11.2 fim PAH feature, the line 
strengths indicate PAH abundances between 10 and lOOx lower 
than the template model abundance, i.e., between 0.06 and 0.6 
% of dust mass. With respect to hydrogen, this translates to 
PAH abundances between 5 x 10~ 9 and 5 x 10~ 8 . Geometry af- 
fects these numbers by at most a factor of a few (see Sect. 15. 5> . 
Thus, the absence of PAH features in the majority of our T 
Tauri sources may be partly caused by a lower PAH abundance 
than found in molecular clouds. However, this does not mean 
that PAHs are absent; even at these lower abundances, PAHs 
can still have an important influence in terms of UV opacities 
and heating rates in the surface layer of the disk (Jonkheid et 
al., submitted). For comparison, the PAH abundance inferred 
for t he HP 141569 tran sitional disk is 0.00035% compared to 
dust (Li & Lunine 2003) or 1.5xl0~ 10 with respect to hydrogen 
Jjonkheid et alJ2006l) . 

5.5. Disk geometry 

Finally, the effect of the disk geometry on the PAH features 
is addressed. The template disk and PAH model is used, again 
with r e ff = 10000 K. The pressure scale height of the disk at 
the outer radius is varied between H p /R out = 0.06 and 0.22, 
while the pressure scale height of the inner rim is kept constant 
at Hp/R[ n - 0.02. This illustrates the effect on the PAH features 
of a flaring disk versus a flatter disk. SEDs for varying values 
of the scale height of the outer disk are shown in Fig. ^2 The 
line flux of the PAH features is found to slowly decrease as the 
disk becomes less flared by lowering the outer disk scale height 
from H/R oM = 0.22 to 0.06. The feature over continuum ratio is 
found to increase slightly with decreasing H/R out for most fea- 
tures, with the exception of the 3.3 and 6.2 yum features, where 
the feature/continuum ratio decreases for H/R out < 0.1. This 
trend is caused by a more rapid decrease of the silicate contin- 
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Fig. 11. Top panel: Model SEDs of a star+disk with PAH, with 
scale height of outer disk H/R out = 0.06, 0.1, 0.14, 0.18 and 
0.22 at a distance of 150 pc. Middle panel: PAH line fluxes 
for the corresponding models. Bottom panel: Feature to con- 
tinuum ratio of the various features. The dotted line indicates 
our observational limits. 



uum flux compared to the PAH peak flux. All the main PAH 
features can be distinguished for the flaring disk, although the 
8.6 fim feature has a low feature to continuum ratio of ~3- 
10%. 

To compare this scenario to our PAH detections, a measure 
for the flaring index of the disk is taken to be the shape of the 
SED , obtained by d ividing the observed vF v at 35 yum over 13 
yum dKessler-Silacci et"a"fll2006l) . This ratio is compared to the 
measured 11.2 /mi line flux. No clear evidence for a correlation 
between the 11.2 fim line strength and a rising/non-rising SED 
is found. Together with the lack of strong changes Fig. ^2 this 
suggests that disk flaring alone is not sufficient to explain the 
PAH non-detections in our sample. 

A second geometry effect, disk thickness, can have a direct 
effect on the PAH feature strength /pah- A geometrically thicker 
disk captures a larger portion of the stellar radiation field than 
a flatter disk, and hence the response in both the dust contin- 
uum and the PAH features is stronger. Observations of spectral 
energy distributions of disks around pre-main-sequence stars 
have indicated that disk geometries vary strongly from source 
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to source dMeeus et all200ll) . and that this factor can therefore 
potentially affect the spectrum. Test model runs show that in- 
creasing the disk geometrical thickness by a factor 2, by in- 
creasing the pressure scale height H p throughout the entire 
disk, introduces an increase in /pah of ~ 1.5-2, while the ra- 
tio of PAH feature over continuum ratio decreases by a factor 
~ 1.4-1.8. The disk would have to be a factor ~ 60-80 thinner 
for the PAH features to have been undetectable for T e ff < 6000 
K (GO), for the template abundance. 

A third geometry effect, inclination, affects both /rah and 
the feature over continuum ratio. Our template models are run 
for face-on disks (i = 5°). When the template model is ob- 
served at angles up to i ~ 60°, /rah remains relatively un- 
changed (< 3%) and the feature over continuum increases by 
about 10%. For an almost edge-on disk (z = 80°), /pah nas de- 
creased by ~ 40%, while the feature over continuum is a factor 
~ 4 larger. Inclination may be an explanation for some of the 
sources with abnormally large F/C ratio, compared with their 
feature strength, e.g. RR Tau. 

A fourth geometry effect which could explain an unusually 
high feature/continuum ratio is the introduction of a large scale 
gap in the inner dust disk, for example as the result of dust co- 
agulation or possibly through clearing by one or more planet(s), 
although the latter would require extreme circumstances. This 
lowers the 3-13 /mi continuum flux while keeping the PAH 
feature strength the same. The motivation for this explanation 
comes from the SED of the small fraction of T Tauri stars in 
our sample that are detected, all of which have a higher 11.2 
/mi feature/continuum ratio than would be consistent with their 
11.2 /mi fluxes within the standard model. The SEDs of these 
4 detected sources show remarkably low emission at 10 /mi 
followed by a steep rise towards 20 /mi, suggestive of a gap 
opening up in the dust disk (Brown et al., in prep.). Examples 
of similar S EDs include those of DM Tau, GM Aur and CoKu 
Tau/4 (e.g. iD'Alessio etafll2005h . for which the continuum 
flux around 11 /mi is also depressed by more than an order 
of magnitude compared with that of a standard T Tauri star in 
Taurus. 

5.6. Model summary 

In summary, all three effects of reducing the stellar temperature 
(and thus the amount of optical and UV radiation), reducing the 
PAH abundance, and decreasing the flaring angle or thickness 
of the disk are shown to decrease the absolute strength of the 
PAH features, as expected. The feature over continuum ratio 
trends are similar for the 6.2, 7.7, 8.6 and 11.2/mi PAH fea- 
tures. That of the 8.6 /mi is the lowest in all models, with typ- 
ical values of < 6-10% even for models which maximize the 
features, which can qualitatively explain the observed lack of 
this feature. 

The 11.2 /mi feature remains one of the strongest features in 
all models explored and this is qualitatively consistent with our 
observations. For the template model parameters which maxi- 
mize the PAH features, PAHs are observable within the mean 
sensitivity limit of our observations for r e jf > 4200 K (K6). 
For a lOx lower PAH abundance, the features are still above 



our sensitivity limit for T e g > 4900 K (K2), but now the fea- 
ture/continuum ratio becomes the limiting factor, with values < 
5% for r e jf < 6000 K (GO). This partly explains the absence of 
any PAH features in the spectra of more than half of the sources 
in the sample considered, which consists mostly of late K and 
M type stars. 

6. Conclusions and future work 

Spitzer-IRS spectra were obtained for a set of 54 pre-main se- 
quence stars with disks, including 38 T Tauri stars and 7 Herbig 
Ae stars. The observations are an order of magnitude more sen- 
sitive than those used in previous surveys of PAHs in disks. We 
detect PAH features in at least 3 T Tauri stars, with an addi- 
tional 14 tentative detections to be confirmed, resulting in a 
lower limit to the PAH detection rate in T Tauri stars of 8 %. 
Spitzer SL observations are needed to confirm the presence of 
the PAH features for sources where we currently only see the 
11.2 ytzm feature in the SH spectrum. All 4 sources that show 
hints for inner holes in their dust disk also a show clear 11.2 //m 
PAH feature. 

The lowest mass source with PAH emission in our sample 
is T Cha with spectral type G8. The derived 1 1 .2 /mi line inten- 
sities are between a few xl0~ 15 and 3.3 xl0~ 16 W irT 2 , which is 
typically an order of magnitude lower than what was observed 
for HAeBe stars with ISO. 

Radiative transfer modeling of disks coupled with PAH 
emission models shows that for stars of late spectral type, the 
11.2 yum feature is expected to be the best tracer of the presence 
of PAHs. The models also show that the 7.7 and 8.6 /mi PAH 
features are most affected by veiling by the continuum due to 
strongly rising silicate emission, resulting in low feature over 
continuum ratios. 

For the small number of T Tauri sources detected as well 
as for the Herbig stars, the measured PAH line fluxes and 
feature/continuum ratios are lower than those found from our 
template disk model which maximizes the PAH emission. 
Variations of the model parameters indicate that the most likely 
explanation is lower PAH abundances by factors of 10-100, 
with geometry affecting these conclusions at the level of a fac- 
tor of a few. The high feature/continuum ratios for the detected 
T Tauri stars are due to their abnormally low continuum at 1 1 
//m caused by the dust holes in their inner disks. 

The template model predictions indicate that the 11.2 //m 
feature strength becomes undetectable at our sensitivity limit 
when T e g < 4200 K (K6). This likely explains the absence of 
PAH features for the majority of sources in our sample which 
have spectral types later than K0 (~ 70%). If a large fraction of 
these sources would have excess UV radiation at the level de- 
tected for some K7 stars, then the absence of PAH features im- 
plies a lower PAH abundance by at least an order of magnitude. 
The same conclusion holds for the 1 1 out of 38 sources with 
spectral types earlier than K6 which are not detected. Geometry 
affects these conclusions by a factor of a few. Thus, the lack of 
PAH detections does not mean that PAHs are absent in these 
disks; even at lower abundances they do need to be considered 
since they can affect disk structure, chemistry and gas heating. 
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In a follow-up paper, the combination of 5-35 j-im spectra 
for the entire sample with IRAC and MIPS photometry will be 
used to improve the uncertainty in the PAH detection rate and 
to compare the PAH sources in terms of their disk SEDs and 
interpret this with the disk modeling. 

Ground-based spectra using instruments like VLT-ISAAC 
and VLT-VISIR on 8-m class telescopes will be able to both 
search for presence of the 3.3/im PAH feature as well as bet- 
ter characterize the shape of the 8.6 and 11.2 //m PAH feature 
through higher spectral resolution. Through higher spatial res- 
olution, they will also permit us to put a stronger constraint on 
the spatial extent of the PAH features from T Tauri disks. 
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Appendix A: Model tests and comparison with 
Habart et al. 2004 

There are various differences in the modeling approach be- 
tween the lHabart et alJ ll2004l) paper and our models, which 
are summarized in Table |A~T1 First, the H04 model is a 1 + 1- 
D model, meaning that the model consists of a series of verti- 
cal 1-D models at different radii, combined to make a full disk 
structure. This method follows the irradiation-angle philos oph y 
used in many oth e r models as well, e . g.lD'Alessio et alJ dl998). 
iBell et alJ i 1 997t) . Dullemo nd et alJ J2002I) . In our model we 
use full 3-D radiative transfer based on an axisymmetric disk 
density structure. Our model therefore also includes emission 
from the dust inner rim, contrary to the approach of H04. Our 
3-D approach allows radiative energy to be exchanged between 
adjacent radii, which is not possible in the 1 + 1D approach. 
Consequently the model SEDs differ somewhat between these 
two types of models. Specifically, for a flaring disk, the radia- 
tion will more readily escape in the polar direction than in the 
radial direction, because of the larger optical depth in the latter 
direction. When the disk is viewed face-on, the SED will be 
boosted in the 3-D approach compared to the 1 + 1-D approach 
by up to ~40%. Moreover, the axisymmetric 3-D models de- 
scribed here allow the treatment of radiative transfer in disk 
geometries that are not flaring, e.g. self-shadowed disks (see 
Section[5]for examples). On the other hand, in our models the 
vertical geometric thickness of the disk is only estimated to be 
roughly consistent with hydrostatic equilibrium, whereas the 



H04 model includes detailed hydrostatic equilibrium. The re- 
sults of these two approaches, however, are nearly identical. 

There are also fundamental differences in the opacities and 
the treatment of dust grains in the radiative transfer. Most im- 
portantly, H04 thermally decouple the carbon grains from the 
silicate grains, while we have both grain types at the same tem- 
perature and mix the opacities into a single opacity. Our as- 
sumption may be more realistic because silicate and carbon 
grains can only have different temperatures at the same loca- 
tion in the disk if such grains are physically disconnected. Only 
a small amount of coagulation is required to get these grains in 
physical contact, forcing them to have the same temperature. 
The result is that the silicate feature strength over the contin- 
uum is much stronger in our case than in the case of H04, as 
can be seen in their Figure 3 as compared to our Fig. [3] If the 
graphite is decoupled from the silicate, the graphite heats up 
and the silicate cools down due to the much lower optical/NIR 
opacity of silicate compared with graphite. The superheated 
graphite emission therefore fills in the spectrum on both sides 
of the silicate feature in the case of thermally decoupled grains. 

For stars with lower effective temperature, radiation at 
longer wavelengths becomes m ore inportant. Our mod els in- 
clude the new PAH opacities of Mattioda et alJ J2005t) at op- 
tical and near-infrared wavelengths, which are higher than 
those adopted by H04. In our model we do not include PAH- 
destruction in the inner disk, in contrast with H04. Test models 
with and without PAH destruction show only small differences, 
however. 

The biggest differences between H04 and our models are 
likely due to the treatment of the stellar radiation. Our models 
assume that the full stellar flux can irradiate the disk, whereas 
H04 assume that the disk goes all the way to the stellar surface 
and occults the lower half of the star. Hence the H04 models 
have only half the illuminating flux as ours. Moreover, they in- 
clude the star as a blackbody emitter with main sequence stel- 
lar parameters, while in our model we use a Kurucz model with 
pre-main sequence values. 

To quantitatively compare the two models, we have plot- 
ted our PAH intensities in the same way as H04, as functions 
of the integrated stellar radiation field in the FUV (6-13.6 eV, 
912-2050 A) wavelength range. The UV flux is usually spec- 
ified rel ative t o the average interstellar radiation field from 
Habing (1968) integrated o ver this wavelength range, which 
is 1.6 x 10" 3 erg cm" 2 s" 1 JTielens & Hollenbacblll985t) . This 
quantity is most often referred to as Go but we refer to it here 
asx to be consistent with H04. H04 assume a single blackbody 
for the stellar radiation field for stars of spectral type B5-G0 
and rescale the flux from the stellar surface R*, to the flux at 
a distance from the star of 150 AU, as a reference point in the 
middle of their disk. 

We have calculate^ for our models as presented in Fig. [9] 
using Kurucz spectra and pre-main sequence parameters. The 
results are presented in Fig. IA.1I Moreover, we have run our 
template model assuming a blackbody stellar spectrum and 
main-sequence stellar parameters, with T e g = 15000, 12300, 
10500, 8500, 7000, 6000 and 5000 K, corresponding to spec- 
tral types B4, B7, AO , A3, F2, F9 and K2, respectively 
llGrav & Corballvll994l) . Fig. lA.2l shows the results for the tern- 
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Fig. A.l. Strength of 1 1 .2 fim PAH feature /pah (scaled to a dis- 
tance of 150 pc) vers us r, the integrated F UV field (6-13. 6eV) 
for a ZAMS star, after lHabart et al.ld2004 their Figure 8). Solid 
line represents our template model, based on Fig. [7J dashed 
and dash-dotted lines represent our models with a lOx and 
lOOx lower PAH abundance respectively. Black '*' symbols 
are for c2d sources. Grey diamonds are ISO upper limits of the 
1 1.2 /mi feature strength for Herbig Ae sources, grey '+' sym- 
bols are ISO detections. c2d sources are labelled as follows. 1: 
TCha; 2: LkHa330; 3: SR21 N; 4: HD 135344; 5: RRTau; 6: 
VV Ser; 7: HD 101412; 8: HD 98922. The dotted grey line in- 
dicates our typical 3<x sensitivity limit of 2.5 x 10~ 16 W rrT 2 for 
sources at d — 150 pc. 



plate PAH abundance (same as H04, dark solid line) and lOx 
lower (dashed line) and includes the model prediction by H04 
(light solid line). The observational data are included in both 
figures, with x computed from the spectral types in the same 
way as described above for the pre-main sequence stars. 

It is seen from Fig. IA. H and Fig. lA.2l that for the same PAH 
abundance, our models predict factors of 4-5 stronger PAH 
features than the H04 models. In fact, our entire SED (contin- 
uum and features) is stronger than in their models. This means 
that our model disk absorbs more stellar radiation and hence 
re-emits more radiation in the infrared, both in PAHs and in 
thermal continuum. A factor of two can be traced back to the 
different assumption for the irradiative flux, whereas another 
~40% is due to the "boosting" effect in the 3-D versus 1 + 1- 
D approach in this near-face-on geometry. The remaining fac- 
tor of ~ 2 is ascribed to the other differences described above 
in terms of disk structure, geometry, treatment of dust grains, 
PAH destruction and PAH opacities. Another difference is that 
our template model assumes a mix of neutral and ionized PAHs 
whereas that of H04 assumes only neutral PAHs. A comparison 
of feature-over-continuum ratios, in which many of these ef- 
fects drop out, would be instructive but H04 do not quote such 
values. 



Fig.A.2. As Fig. IA. II but assuming a blackbody stellar spec- 
trum and main-sequence stellar parameters. Dark solid line rep- 
resents our template model, dashed line our model with a lOx 
lower PAH abundance. Light solid line represents the H04 disk 
model results for the 11.2 fim feature. 

Comparison of Fig. lA.ll and Fig. IA. 21 shows that the 1 1.2 
fim fluxes are lower if blackbody radiation is assumed with 
main sequence stellar parameters, especially for T Tauri stars 
(later than F8). Both our blackbody models with lOx lower 
abundance as well as the H04 models fit most of the ISO and 
Spitzer detections for HAeBe stars. Both models predict that 
for T Tauri stars of spectral type G5 and later, /pah to falls 
below our Spitzer detection limit for d — 150 pc. However, 
the H04 models underpredict the detected 11.2 fim strengths of 
some of the lowest luminosity sources in our sample, whereas 
our models with standard PAH abundance come much closer. 
In general, the PAH emission from T Tauri stars, when de- 
tected, appear better modeled with Kurucz spectra for pre-main 
sequence stars than with blackbody spectra for main sequence 
parameters. Overall, the main conclusions from our paper are 
not affected qualitatively if the H04 models with blackbodies 
were used to interpret the data. Quantitatively, the threshold for 
detection shifts to earlier spectral types (G5 versus K5) and the 
inferred PAH abundances are factors of 4-5 higher. 

Finally, we note that a multitude of tests of our code were 
performed in addition t o the tests that the RAD MC code already 
has gone through (e.g. Pas cucci et alJl2004l) . We tested that 
our model conserves luminosity; we compared the SED with- 
out PAHs to the SED prod uced by the C hiang & Goldreich- 
type model described in Dullemond et all (I2OO 1 ): we estimated 
what the SED strength should be on the basis of the disk geo- 
metric thickness; and we performed optically thin models con- 
sisting purely of PA H grains and compared to the models of 
iLi & Drainel J200lh . All these tests confirmed the validity of 
the models. 
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Table A.l. List of differences between the model presented here and that of Habart et al 1 d2004h . 



IVToHpI fisnppt 


This nanpr 


H04 


Radiative transfer 


Axisymmetric 3-D 


1 + 1-D 


Disk structure 


Inner rim 




P(Z) 


Parametrized 


Self-consistent 


Stellar spectrum 


Kurucz 


BB 


Stellar parameters 


PMS 


MS 


Stellar flux 


100% 


50% 


PAH evaporation 


1/2 


yes 


T -decoupling 


no 


yes 


Optical-IR PAH opacities 


Li &Draine (20011. 


Li & Draine (2001). 




Mattioda et al. (2005) 


Desert et al. (1990) 
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Appendix B: Online Material 

The Online Material contains Figs. IB. ll and IB. 21 presenting 
blow-ups of the 11.2 and 12.8 pm PAH features. 
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Fig. B.l. Left: Blow-up of the Spitzer IRS high resolution spectra around the 11.2 fim PAH feature. A simple fit of the continuum 
flux below the PAH feature is plotted with a dotted line. Right: Continuum subtracted spectra, normalised to the peak flux of the 
PAH feature. 
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Fig. B.2. Left: Blow-up of the Spitzer IRS high resolution spectra around the 12.8 fim PAH feature. A simple fit of the continuum 
flux below the PAH feature is plotted with a dotted line. Right: Continuum subtracted spectra, normalised to the peak flux of the 
PAH feature. 



